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The nonadiabatic correction to spin transfer torque arising from fast-varying spin texture 
is calculated treating the conductions electron fully quantum mechanically. The torque is 
nonlocal in space, and is shown to be equivalent to a force (due to momentum transfer) 
acting on the center of mass of the texture. Another kind of force exists in the adiabatic 
regime, and it is identified to be of topological origin. These forces are shown to be the 
counter-reactions of electric transport properties, resistivity and the Hall effect. 
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1. Introduction 

Recently, torques on local spin (magnetization) induced by spin-polarized electric current 
and spin current are of special interest from the viewpoint of magnetization switching by 
current in wires 1-3 and nanopillars. 4 ' 5 Berger 2,5 and Slonczewski 4 pointed out that the ex- 
change of the spin angular momentum between the local spin and current-carrying conduction 
electron induces a torque called spin-transfer torque, and that this is a dominant mechanism 
of magnetization flip when the magnetization is slowy varying (i.e. in the adiabatic regime). 
In the case of a domain wall, another effect arising from the current, that is, the effect due 
to elecron reflection or nonadiabaticity, has been discussed by Berger. 1 This force was later 
derived microscopically by estimating the nonadiabatic correction of the torque due to ex- 
change coupling. 6 Recently, a new torque term in the Landau-Lifshitz(-Giibert) equation was 
proposed, which was shown to arise from spin relaxation of the conduction electron. 7,8 This 
term, sometimes called the /3-term, was later derived microscopically 9, n ' 12 and also from other 
mechanisms. 13 This term is also sometimes called nonadiabatic torque, in that the term arises 
from the deviation of the electron spin from perfect adiabaticity as a result of spin relaxation. 
(This effect of spin relaxation is not the nonadiabaticity in the strict sense, so we will call it 
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" nonadiabaticity due to spin relaxation" to avoid confusion.) 

Quite recently, Waintal and Viret 14 and Stiles and coworkers 15, 16 studied the spatial dis- 
tribution of the current-induced torque around a domain wall by solving the Schrodinger 
equation and found a nonlocal oscillation torque. This torque is due to the nonadiabaticity 
arising from the finite domain wall width, or in other words, from the fast-varying compo- 
nent of spin texture. The oscillation is of period ~ /cf™ 1 (fe.fr is the Fermi wavelength) and 
is of quantum origin similar to Ruderman-Kittel-Kasuya-Yosida(RKKY) oscillation. The os- 
cillating torque is due to oscillating spin accumulation of the electrons as argued in ref. 17 
Quite recently, nonlocal oscillating torque around domain wall was numerically studied tak- 
ing account of strong spin-orbit interaction on the basis of the Kohn-Luttinger model (i.e., in 
magnetic semiconductors). 18 It was shown that the oscillating torque is asymmetric around 
the domain wall and that this feature results in high wall velocity. 

These results can be summarized as a total torque acting on local spin, S, given by 

T = ~^s Us ' v)s ~^i [Sx {js • v)s] + Tnl ' (1) 

where the first term represents the spin transfer torque, the second is the (3 term due to spin 
relaxation, and the last term denotes the nonlocal torque argued in refs. 14-16 (but explicit 
expression was not obtained there). We note here that the non-adiabaticity has two different 
origins, spin relaxation and fast-varying local spin texture, and that these two lead to to- 
tally different torques ((3 term and r n i). Spin transfer torque was studied further within the 
Boltzmann approach in ref. 19 So far, most numerical simulations on magnetization dynamics 8 
have modeled the effect of current solely by local torque terms and nonlocal torque r n i has 
been neglected. Such simulations are thus valid only in the slowly varying limit. Although 
the nonlocal term is believed to be small in most cases, this term could still be important in 
reality, since, like the (3 term it plays a qualitatively different role as spin transfer torque, and 
f3 term is also quite small. In particular, the nonlocal term would be important in vortex-wall 
structures, where the spins vary rapidly inside the core. As we will demonstrate, two nonadi- 
abaticities, (3 and r n i, both contribute to a total force acting on spin texture. The dynamics 
of a rigid domain wall under the total torque given by eq. (1) was studied in ref. 20 

One of the aims of this study is to derive an explicit expression of nonlocal torque, r n i, 
within the lowest order correction to the slow variation, treating the conduction electron 
fully quantum mechanically. As we will show, quantum treatment is essential for discussing 
nonadiabatic torque and force. In our calculation here, we assume Ar/h » 1, where A and r 
are exchange splitting and (elastic) lifetime of the conduction electron, respectively, and takes 
account of the nonadiabaticity represented as higher order terms in l/(kp\), A being the scale 
of the spatial variation of spin texture. Ar/h » 1 is one definition of adiabaticity 21 ' 22 where 
the gauge field can be treated perturbatively. 
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Second, we will show that this complicated nonlocal torque represents the effect of electron 
scattering by the spin texture and that the nonlocal torque is summed up simply to a total 
force acting on the texture. Third, this force is shown to be exactly proportional to resistivity 
due to spin structure and hence to the reflection probability. Finally, we will show that there 
exists another type of force which remains finite in the adiabatic limit. This force, adiabatic 
force or topological force, is proportional to the topological vortex number of the texture, as 
was pointed out in ref. 10 ' 23 It is shown to be proportional to the Hall resistivity caused by 
spin chirality or vortex and is thus identified with a back-reaction of the Hall effect due to spin 
chirality or the spin Berry phase. 24-28 This force acts on a vortex or vortex wall perpendicular 
to the current. The effect of this topological force is included in the conventional spin transfer 
torque term, in contrast to the nonlocal torque. 

This paper is an extension of the preceeding paper. 6 While torque and force acting on 
a domain wall were exclusively studied there, those for general spin structures are studied 
in the present paper. Besides, in ref., 6 spin transfer torque was approximated by dominant 
contribution (i.e., in the adiabatic limit) and non-adiabaticity was included only as a force. 
This inconsistency is removed in the present paper, calculating both to the same order of non- 
adiabaticity. Another aim of the present paper is to calculate the effect of current as the linear 
response to the applied electric field. This was not done correctly in ref., 6 where the effect of 
current was assumed to be represented by the distribution function of current-carrying state. 

2. Model and Method 

The Hamiltonian we consider is the standard exchange interaction one given by 

H e = ^ ^ ^kcr^ kcr c ka "i" -ffimp; (2) 
kcr 

where electron operators are denoted by Cfc CT and c\ a , = ^ — £f, cr = ± denotes spin, and 
#imp represents the elastic scattering by impurities (without spin flip). Conduction electrons 
are coupled to local spin, S, via the exchange interaction 

H cx = ~ J d 3 xS(x) • J x ac Xl (3) 

where <x is the Pauli spin-matrix vector, and A is the exchange coupling constant. The La- 
grangian of the conduction electron is defined as L c = h f cPxic'c — H c . 

All the torque acting on the local spin caused by the electron is thus obtained from the 
exchange term as 

T{X) = J§§) X S{X) = ~P {X) X S{X) > (4) 

where 

s(x) = (c ] x ac x ^ (5) 
is the spin polarization of the conduction electron evaluated in the presence of spin texture, 



3/25 



J. Phys. Soc. Jpn. Full Paper 

S. Our task is merely to carefully estimate this spin density. 

To incorporate the spatially varying spin texture, we consider the slowly varying case, 
called the adiabatic case. In this case, gauge transformation in spin space is useful. We consider 
only the lowest-order contribution from the gauge field, but will go beyond the adiabatic limit 
by taking account of the nonadiabaticity represented by the finite momentum transfer {q 
defined below) between the electron and gauge field. Strictly, this expansion is justified in the 
limit of Ar/^> l. 21 ' 22 

Gauge transformation, defined by a 2 x 2 matrix, U(x,t), relates the electron operator c 
to a new operator a as 

c(x, t) = U(x,t)a(x, t), (6) 

where U is defined as 

U(x,t) = m(x,t) ■ a, (7) 

m being a real three-component unit vector. The matrix satisfies U 2 = 1, or U^ 1 = U. The 
aim of our gauge transform is to let this spin to be along the z-axis, i.e., ^U^(S ■ a)U = 
2m(m ■ a) — a = a ■ e z . This is satisfied if we choose 

( , 6 , 8\ 
m = I an - cos <p, sin - sin <p, cos - I . (8) 

We see that a gauge field appears from the derivative terms, as seen in 

<V = U(d^ + U- l d^U)a = U (<9 M + iAJa, (9) 

where a gauge field is defined as 

A„ = -iU~%U. (10) 

In terms of m, is written as 

A fl = (mx d^m) ■ a = A*a a , (11) 

where the summation over a = x, y, z is suppressed. 

The electron part of the Lagrangian is written in terms of the a-electron as 



ifuJa - e ka (J ka a ka - ft£ (j^ (k + | ) • A°(-q)j a j k+q a a a k 



2m 

QP 



(12) 



where 



£fc CT = ^J— - £F - crA (13) 
2m 
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is the electron energy in a uniform field along the z-axis and 

J„(fc)= (14) 

for fi = x,y, z, t. In this paper, we assume slowly varying in time and neglect the t- 
dependence. 

The electron spin is modified by the gauge transformation and is written using the spin 
density in the gauge-transformed frame, s, as 

s(x,t) = (c'trc) = 2m(m ■ s) — s, (15) 



where 

s(x,t) = (attra^ . (16) 
The original electron spin is written as 

s = -s\\e + s z n - s±e ( j ) , (17) 

where sm and s± are defined as 



2 



e+ p s = cos q>s x + sin <f>s y 



s± = - ^2(^f)ie Tl< ^s ± = - sm 4>s x + cos 4>s y , (18) 



± 

with = s x ± is y . The three unit vectors are n = S/S = (sin 6> cos c/>, sin sin 0, cos 6), 

eg = (cos 9 cos 4>, cos 6* sin 0, — sin 9) (19) 
= (— sinc/>, cos c/>, 0) . (20) 

3. Spin Density in Linear Response 

We will estimate spin density s(x, t) to the lowest order in the gauge field. Electron spin 
is modified by spin texture even in the absence of current. (This effect was not considered in 
ref. 6 ). This equilibrium contribution is given by (Fig. 1) 




Fig. 1. Elecron spin density without current at the lowest order in gauge field (represented by wavy 
lines). 
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5±(0) 



kq 

X [^-f, T (^)5fc + £, ± H+5fc_£, T (^K + £ ± H] 

= y E e /fc g +f ~ h ~ h l i » ^ 21 ) 

fcq fc+2",± fc— t 

where g ka (u) = 1 +J _ is a free retarded Green function and g < denotes the Keldysh 
Green function. 

Let us next calculate a current-driven part. We estimate the response to a static applied 
electric field to the linear order. The interaction between the electron and applied field is given 
as 

Hem = J d 3 xA EM -j, (22) 

where 

A EM = (23) 

is the electromagnetic gauge field, E is the applied electric field assumed to be spatially 
homogeneous, and Q,q is the frequency of the applied field, which is chosen to be Qq — ► at 
the last stage of the calculation. We note that this linear response calculation takes care of the 
current correctly, while treatment described in ref. 6 lacks consisitency in that a generalized 
Fermi distribution function (/*.) in the presence of current flow was assumed there. The total 
electric current is given in the gauge-transformed space as 

/,o eh i , j. *± \ eft v-^ + 
dx--(J V z c)-—A EM ^clc k 

k 

= ^ E ka i ak - ^ E A >i+ q aa<1 k - ^ E A ™ a W- (24) 

k kqa k 

Spin density at the first order in both of gauge fields, Aem and A^, is shown in Fig. 2. 
The first two processes (denoted by (A)) are obtained as 

Bi 



S ± (x,t)^ = EE^[( fc + |)j5 fc _ fjTi _^5 fc+fi± ^£| i ff fc+fi±iW+ £| 1 ] < 

u)kq i 

+ ( k " I), [9 *-i ,=F,«+^*+§ ,±,"+^ ]< ] ' (25) 



where 

R = 

and the last contribution (B) arising from the modification of the vertex is given by 

— n r>„ n o„ 

^0 



Bt = ^ e -*(*-no«) £ J M (fc)A± (g), (26) 



uikq 
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Fig. 2. Diagrams representing the current-driven part of the electron spin density at the linear order 
in both gauge field (represented by wavy lines) and applied electric field (dotted lines). Upper two 
processes are self-energy corrections and the lower process is from the correction of the current 
vertex. 



where 



Bn = 



eEj 
mV 



e- i{q - x -^Af(q). 



(28) 



We will derive the leading contribution in the clean limit, r — > oo. The calculation is lengthy, 
so we present simply the result here. Details are presented in the appendix. The result is 
obtained as (s±W = 5±( A ) + s ± ( B )): 



5±W, 



x,t) 



j_yy^ e - iq . x y A ± {)T 



2vr ^ ^ 2mV 

kq i 



(-(*+!)) 



Q\\ 9 k ,±-9 k ,± 



_ _ k .j ( k , 1) 9 k , T -9 k , T 

e k+g -e k ±2A 1 »\ 2Je k+q -e k T2A 



(ai, T Jn (k + |) 5{e k+q - efc =F 2A) - ^ ± J M (- (k + |)) <5(e fc+g - e fc ± 2A) 



(29) 



In terms of sy and sj_, the equillibrium contribution (eq. (21)) is written as 
s m (x,t) = ^-Y. e ~ l<l ' X [%>(Q,*)Hi(q) + A a (q,x)H 2 (q)] 



_(0) 



= IE 6 " 19 " ^o(<7,*)#i(<?) -^(<jf,aj)ff 2 (g)], 



(30) 



where the correlation functions are given by 

~v 

H 2 (q) 



H M = ^E p fk± 



k± ±2A + ^! 



fc± 



2k-q + (t 

2m 



(31) 
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and A s ^(q,x) = | E± e^^A* and A*(g,x) = \ E±(T«)e Ti0(a:) Aj (they depend on both 
q and x). These can be written as Af(g, a;) = e S( p(x) ■ A )1 (q) and A"(g, x) = e^x) ■ A^(q), 
where e S( k = (cos 4>, sin eft, 0). Integration over fc is carried out to obtain 



TT ( \ _ 3 1 1 X ^ 



^(g 2 -i)(? 2 -C 2 )in 



(g + l)(g±C) 



(g-i)(g=FC) 



±(i±C)(C±g 2 ) 



(32) 



where g = |<jf|/(2fc^), s = ^rC(3 + C 2 ) is the electron spin density, k F = r\(kp + + kp- 



C = fc^ + + ^ i an d Hi is normalized to be H\ = 1 + 0(q 2 ). We also obtain 



H 2 (q) 



k F 3 ( l 

16"7T g 



i+C 



dkk6{kq - \q 2 ± C)|) = k F 3 H 2 (q), 



± -/i-c 



(33) 



where the integral is to be taken only in the regime k > ||(<T i C) I- Correlation functions, -ffi 
and H2, are plotted in g-space and in real space in Fig. 3. 




Fig. 3. (Color online) Plot of correlation functions, H\{q) and ^2(9), and their Fourier transforms, 
Hi(x) and H2{x), describing nonlocal component of spin density (and torque) in the absence of 
current. H\ has a finite adiabatic component (q = 0, i.e., local component), while H2 does not. 
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Similarly, current contributions are calculated as 

sWW(x,t) = ^ e -^[(^-i s )xi(q) + (^-i a )X2(Q)] 

~s ±W (x,t) = ^Te-^ x [(E.A*) Xl ( q )-(E.A s ) X2 (q)], (34) 

where 

X2lq) = £^T,{ ± ^)^i k + W s(Ck+ '- Ck±2msii - g% ^ (35> 

fc± 

The summation over k in eq. (35) is carried out using y ^ fc = J dev(e) d c ° s ® k , where 
6k is the angle of k measured from the q direction. We carry out energy integration first, 
assuming that the dominant pole arises from the green function and the residue contribution 
from 1 i 0A is neglected. After some calculation, we otain 

6T 

Xi (<?) = — (n+-n_)xi(g) 



3 + C 2 ± 



2 i /• z2 



1 H ; in 

2(1 ± Qq 



;i + 9)(9±o 



(1-?)(9=FC) 



(36) 



where n± = -^fr is the spin-resolved electron density. \ ls normalized to be x = 1 + 0(q 2 ) 



and thus Exi(q) = -k?'sXi(<7) {js = E(n + — n^)e 2 T/m is spin current vector). We similarly 



obtain 



where 



and 



m 2 A 2 (Jo M<?) CT o - , x , QV x 

X2(9) = "I^T— = T X2((?) ' (37) 



~ / \ c 2 Mg) , Qe x 

X2(9) = - i l + 3C 2 g ~ . (38) 



M«) = < . (39) 

[ otherwise 

represents the regime of Stoner excitation, do = e 2 nr/m is the Boltzmann conductivity, and 
n = n + +n- is the total electron density. As is obvious, for small q, \2 = 0. In the real space 
representation, 

SU1) W = ^E /^[Jlxi(^-^0(^N-A l (^0)+/X2^-^(^N-A,(^))] 

i 

(40) 
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= ~^Y^ I ' d3x ' i^Xi( x - x'){e^(x) ■ Ai(x')) - fx 2 (x - x')(e s $(x) ■ Ai(x'))] . 

(41) 

Correlation functions, xi an d X2 5 are ploted in g-space and real space in Fig. 4. Note that xi- 




2k p x 



2k x 

F 



Fig. 4. (Color online) Plot of correlation functions Xi(q), X2(o) and their Fourier transforms, Xi(x) 
and \2(x), describing nonlocal components of spin density and torque in the presence of current. 
Xx has a finite adiabatic component (q = i.e., local component) while X2 does not. 



and X2-terms are proportional to spin current and charge current, respectively, only in the 
adiabatic limit (q = 0), but are not necessarily so when nonadiabaticity sets in, since Xi(q) 
and X2(q) can depend on polarization £ in a complicated manner. 

4. Torque in Landau-Lifshitz-Gilbert Equation 

The electron spin density induces a torque on local spin at x, given by 

-—(S x s) = A(-s ± e e + s^) (42) 



T[X 
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Namely, 8- and (^-components of torque are respectively given by 

Tg = -A S ± 

r = AS,,. (43) 

Thus, if one solves the Landau-Lifshitz(-Gilbert) equation of local spin taking account of 
eq. (42), one can include all the effects arising from conduction electrons. From the above 
calculation, we now have full spin density as s u = s u ^ + s u ^ [y =||, _L). The torque acting 
on a spin at x is then obtained as 

T = T (o) +T (i) j (44) 

where is an equillibrium part, 

r^(x) = J d 3 x [H^x - x'^e^x^e^x) ■ A (x')) - e e {x){e^{x) ■ A (x'))) 

+H 2 (x - x'^e^x^e^x) ■ A (x')) + e e (x)(e s4> (x) ■ A (x')))] . (45) 
and the second part is due to current: 
T W {x) = J2 J j j d 3 x [xi(x - x')[ee{x){ e(t> {x) ■ A^x')) - e^x){e S(j> {x) ■ Mx'))]] 

i 

+ J2 7 / d * x> & 2 ( x - x')[e e (x)(e s4> (x) ■ Ai{x')) + e^x^e^x) ■ Mx')]] . 

i 

(46) 

In the adiabatic limit, H\{x — x') — ► s5(x — x') and x\(x — x') — > S(x — x') become local, 
and H2 and \2 vanish. Thus, total torque reduces to a local one, 



r^x) = 



e e (x) x eJx) x 



sA (x) + *yJ—Ai{x 



If a A 

2 1 sdt + T Js ' v ' n ' 



(47) 



where we use (e sif> (x) ■ A^(x)) = (e e (x) ■ A^x)) = A\ and (e^(as) • A^x)) = A^. When 
nonadiabaticity sets in, we see that (1). the torque becomes nonlocal, (2). torques in 6- and 
0-directions become mixed, as seen by comparing \i an d X2- The second role of mixing 
torque components is very important, since it means that the nonadiabaticity mixes spin- 
transfer torque (r 2 ) and force (F). The oscillating torque found here is consistent with previous 
observations. 14 ' 15 The oscillation is of the order of 2fcp, as seen in Fig. 4. Thus, this nonlocal 
torque is of quantum origin. The classical estimate of conduction electron spin or torque (such 
as on the basis of the Landau-Lifshitz(-Gilbert) equation) thus fails when nonadiabaticity is 
to taken into account. 
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Let us denote the torque above as 

t = r ad + r n i, (48) 

where r n \(= r — T a d) is purely nonlocal. We here take account of torque arising from spin 
relaxation as considered in ref., 11 

3 

where a s f and /? s f are proportional to the spin relaxation rate (see eqs. (82) and (83) below). 
(The spin relaxation process also produces nonlocal torque, but this contribution has so far 
not been addressed in metals, since it would be a nominal small correction. See also the 
discussion at the end of this section.) The total torque due to electrons so far known is 
therefore summarized as 

Ttot = T ad + T n l + T s f , (50) 

and the effective equation of motion of local spin is given by 

d t S=^SxS + ^B cS xS + r tot , (51) 

where «o and B e ff are Gilbert damping and the effective field, both from nonelectron sources. 
Using eq. (47), we can rewrite eq. (51) as 



3 

n = |S x S + £s&t(S x (j s • V)S) + r nI + ^B efr x S, (52) 



where a = a>o + a s f- We see here in the first term that the magnetization is now S + 1 , i.e., the 
total spin, comprising local spin and conduction electron spin. The importance of using total 
spin instead of local spin in current-driven domain wall motion was stressed recently, 13 but 
extra magnetization from conduction electrons is naturally taken account of in the present 
formalism (by a time component of the gauge field). We stress here that this effective Landau- 
Lifshitz-Gilbert equation (with (3 and nonlocal torques) was mathematically derived. Thus 
Gilbert damping, which was recently argued to be incorrect in the presence of current, 29 is 
justified. As we see here, a and (5 are finite and does not necessarily equal, 11 in contrast to 
the observation in ref. 9 Very recent study 12 also supports non equal a and (3.. 

Nonlocal oscillating torque around the domain wall in magnetic semiconductors was nu- 
merically studied quite recently. 18 It was shown that the torque is asymmetric around the 
domain wall because of strong spin-orbit interaction and that this interesting feature results 
in fast domain wall velocity. Such an effect would be described by the spin-orbit correction to 
r n i, but this has not been addressed so far in cases of metals (where spin-orbit interaction is 
caused by random impurities). 
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Let us derive the expression for the force acting on the spin structure. The force is defined 
as the energy change when the structure is shifted. Considering a displacement x — > x + X, 
the force due to electrons is 1,6 

Using ^VS = egVO + e^sinOVcj) and n x eg = — e^, n x = eg, we obtain 

F = 2A J d 3 x (s||A ± - s ± A^ . (54) 

(This expression is equivalent to that described in ref. 23 ) From eqs. (18) and (29) (noting 
g r k(7 ~ — iir5(eka) on the fc-integral) we obtain 

± kq ij 

-ir« (fe + |) S(i k+ ,-e h ±2A)(^Af{q)A : f(-q)+Aj(q)A ± (-q)") . (55) 

The second term vanishes in the adiabatic limit (q = 0) while the first term remains finite. 
We now demonstrate that the second term represents the reflection of the electron due to the 
spin texture. In fact, we find, by using e k + q — e& = ■ (k + f ), that the second term (we 
call F( na )) is 

^(na) = 47TC ^y A2 S(e k± )S(e k+q -e k ± 2A)Af( q )Aj(- q ). (56) 

± kq i 

We can show that this force is proportional to the resitivity due to spin structure, which is 
obained, on the basis of Mori formula, as (in the case of current along the z-direction) 

Ps = Yl \ Aa z(<l)\ 2 ^k +q ,-a ~ ek,*)S{e k , a ). (57) 

kqa 

In fact, 

F (na) = l^L p ^ V = eNcPsj} (58) 
m 

where N e = nV is the total electron number. This is the result in ref. 6 extended to a general 
spin structure. The derivation of resistivity using the Mori formula is shown in §6. 

If reflection is included in F^ na \ what then is the origin of the first term of eq. (55)? This 
becomes clear if we consider the adiabatic limit. In fact, in this limit, the first term reduces 
to 

± kq i 
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• Jsi 



I 
e 



J d*x (A+(x)Aj(x) - AT(x)A+{x)) . (59) 



By using 



and 



we see that 



where 



A+{x)A j {x) - A-(x)A+(x)) = -dsme(di6d j( j> - fyflfl^), (60) 



diS xdjS = nS 2 sin 9^98^ - d^d^), (61) 

Ff = ^ E v / d ' xS • ^ S x 9 > s) = 2n E t 1 '""- (62) 



% = ^3 y ^ 3 xS • (5,S x djS) (63) 
is a vortex number defined in three dimensions. In the case of a thin system, this reduces to 

$ij = n v d, (64) 
where d is the thickness of the system, and 

n v = ^3 J d 2 xS ■ (d t S x djS), (65) 

is a topological number. This force is, in fact, a back reaction of the Hall effect due to spin 
chirality, 24 " 26 and was derived first phenomenologically by Berger, 30 then rigorously in ref., 23 
and also assuming a vortex structure in ref. 31 Note that the adiabatic force is included in (adi- 
abatic) spin transfer torque, while nonadiabatic force is not. The two forces are schematically 
described in Fig. 5. 




Fig. 5. (Color online) Schematic summary of three forces acting on spin structures. Electron reflection 
pushes the structure along the current flow (F na ), and Hall effect due to spin chirality results in 
a force in the perpendicular direction (F ad ). Spin relaxation results in a force along the direction 
of current (F 13 ). 
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6. Resistivity and Hall Resistivity Calculated Using Mori Formula 

The resistivity of the domain wall was calculated by Cabrera and Falicov 32 and later 
by others. 33-36 Here, we derive the resistivity (both diagonal and off-diagonal (i.e., Hall) 
components) due to the general spin texture on the basis of the linear response theory (Mori 
formula 37 ' 38 ) 36 in the clean limit of r — ► oo. The Mori formula relates the resistivity psy (ij 
being spatial directions) to the correlation of random forces in the weak scattering case as 

f e 2 n\~ 2 1 

PS ^ = {-^) S^fa^kv^-xuM- (66) 

Here, Xj jX^e) = ~~ (h/PV) < Ji(i^e)Jj(~^e) > with Jj = dJi/dt = ^[H,Ji], where H is 
the total Hamiltonian and Jj = ^^2 k kid^ k c k ^ s ^ ne total current. The correlation function 
Xjj{fuv) in eq. (66) denotes an analytical continuation of the correlation function calculated 
for imaginary frequency, i.e., Xjj(^) = Xjji^i ~^ ^ + *0). 




Fig. 6. Diagram representing the resistivity in Mori formula at the lowest order in the gauge field. 



The nonconservation of the current (i.e., finite j) arises from the scattering by the spin 
texture. In fact, eq. (12) leads to 



J, = 4 (i)£ 

k,q 



-2A^A. a a\ +q a a a k 



h 2 



+ — (Af((2k + q)-q)- qi(A a ■ (2k + q)) - h qi A%) a) ' k+q a a a k 



(67) 



where we neglect higher order terms in A. We consider a static spin texture (i.e., A$ = 0), 
and then 



Ji = ~ i2 (^) A £ £ A a a ] k+q a a a k . 



(68) 



k,q <y 



The Fourier transform in imaginary time r, defined as jj(iu^) = f^e luJl1 jj(r), where uJe = 
^ is a bosonic thermal frequency, is given as 

Ji(iuj e ) = -VI [^j A ^ ^2 A i C T<l] ^ k+q,n+e<7aak,n, (69) 



k,q,n cr 



where a k<n = e luJnT a k (r)dT and uj n represents the fermionic frequency u n = n ( 2r J 3 ^ , 
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The thermal correlation function is then obtained as 

/2eA\ 2 1 



(2c A \ 1 1 



fcgcr 

where the imaginary time Green function is defined as 

1 



lLO r , 



The summation over uj n is carried out to obtain 



V 7 fcgcr 



Thus the resistivity is obtained as 

4A 2 1 



psa = lim -^-^t — 

r tJ u^o e 2 n 2 V 



E Im 



kqa 



AUq)A-°(-q] 



1 f(ek+ g , 



/(<*,, 



w e fc+ g_ CT - e fe)(T -u - iO 



Then, choosing the current direction as z, resistivity is rewritten as 

4vrA 2 1 



PS = PSzz 



e 2 n 2 V 



E \ A<J z(<l)\ 2 ^k+ q ,-a ~ e k>a )5{e kta ). 



(70) 



(71) 



(72) 



(73) 



(74) 



kqa 



Thus, the reflection force of eq. (56) is proportional to the resistivity due to spin, as we 
explained for eq. (58). The fc-summation can be carried out easily as 



2k Fa q 



0st(q), 



(75) 



where q = \q\. 

Resistance, Rs 



hps-, due to spin texture is related to the reflection probability R, 



according to four-terminal Landauer-Biittiker formula, as i?s = fry^R) an d hence eq. (58) 
indeed relates the force to the reflection of the electron. Equation (58) can also be written, 
using the density of states N(0) = (neglecting spin splitting), as 



JT- = i.w»(0)»„»£* 



(76) 



where Vs = Rsl is a voltage drop by spin texture and / = Aj is current (A being the cross 
section). This equation clearly indicates that the force is due to the momentum transfer (2/cpft 
per electron, with frequency of ^^) multiplied by the number of electrons that contribute to 
the resistance (|eVsiV(0)). 

In contrast to the diagonal component of resistivity, which vanishes in the adiabatic limit, 
the off-diagonal (i.e., Hall) component remains finite in this limit. In fact, we find 



^ ad 
PSij 



4A 2 



T7E Im 



kqa 



A°{q)A-°{-q) 



1 /(<*,- 



/(efc,«r) 



2a A - u 



he 2 



^E Im K(^7 CT (-^)]E/(^) 

q ka 
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n + — n_ 



/2tt c ! z P 
#x». W 5xa,a) = w .,. (77) 

(Correctly, we have retained the antisymmetric component, \{psij — PSji), m the second line.) 

This is the Hall effect caused by the spin chirality, or spin Berry phase, 24 demonstrated in the 

slowly varying case. 25,26 Comparing eqs. (62) and (77), we see that the force in the adiabatic 

limit is exactly due to the Hall effect from spin chirality: 

2e 2 ra 2 X - . 
J n + — n_ ^-r ' J 



= e 2 nJ2j iPSij . (78) 

i 

7. Spin Relaxation 

Let us briefly look into the role of spin relaxation discussed in ref. 11 The spin flip scattering 
was introduced there as 

H a = u s j d 3 x Si5{x - Ri){Jac) x , (79) 

where Si represents an impurity spin at site R4. (The spin-orbit interaction leads to essentially 
the same results as the spin flip case here.) The quenched average for the impurity spin was 
taken to be 

a i S 2 (a,3 = x,y) , , 

SfS = Si jSa ,x\ J_ P (80) 

Kohno et al. 11 considered the adiabatic limit, and obtained the torque due to spin flip scat- 
tering as 

r sI = ySxS| ^3 si (S x (j s ■ V)S), (81) 
where coefficients are given by 

' 2& z v + v- +Sl(v 2 + + » 2 )] . (82) 



a s { = Tm s u 2 



{Sl + Sl)v + + Usl -S*) V - 



(83) 



where v± = m 2 ^ ± and n s denotes the concentration of magnetic impurity. The spin relaxation 
process also produces nonlocal torque, but this contribution has not been addressed so far, 
since it is a nominal small correction in the present context. 

The torque, r s f , is due to the spin polarization caused by spin relaxation (r s f — — i^S xs s f), 

Ssi = ~ ( asf ^ + is Psi {is ' v)s ) ■ (84) 

This spin density gives rise to a force, 

Ff = ^J ^V t S • (a*S + ^PAJs • V)s) • (85) 
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The force due to the (3 term is thus given as 7 ' 8,23 

^ = i 7/?sf js ' (86) 
where 7 = J d 3 x(ViS) 2 , with i \\ j s . If we consider a quasi-two-dimensional system, the total 
force due to current, F = F( na ) + _p( ad ) + F@, is summarized as 

F=j (eN ePs + P-y^Pd) + P^j x g, (87) 

where P = j s /j is the polarization of current and g = ^-<b xy e z is the gyrovector of vor- 
ticity. 31 ' 40 We see that the two forms of nonadiabaticity, due to fast spin texture and spin 
relaxation, both result in a force in the same direction, that of current. 

8. Motion of a Vortex 



Let us consider specifically the case of a single vortex in a film in a xy-plane. As was 
explained by Shibata et al. 31 the Lagrangian of a single vortex in terms of its position, 
X = (x,y), is given as 

L v = l -g ■ (X x X) - U(X), (88) 

where g = 47r {p & xy , and U denotes all the potential energy, including the effect of current. 
(Note that they considered the effect of current in the adiabatic limit separately from U as a 
spin-torque term.) The equation of motion is thus given by 31 

X = ie z x (F-aX), (89) 

where F = -^U(X) and a = a^f- d j d 2 x[(V9) 2 + sin 2 #(V0) 2 ]. As seen, vortex has a velocity 
perpendicular to the applied force, X _L F, which might sound strange but is an interesting 
feature of a topological object. The force due to current is given by eq. (87), and so eq. (89) 
reduces to 

where / = 8n g2$ — (2Se 2 N c ps + jP(3 s f). The perpendicular force on the vortex of the adiabatic 
origin (j x g in eq. (87)) is thus simply a spin torque, which induces the vortex to flow in 
the current direction. This is easily understood since spin torque and the time-derivative term 
of the Landau-Lifshitz equation are combined into a Lagrange derivative along current flow, 
dt — > (pt — Js • V) , which indicates that spin transfer torque drives any spin structure 
along the spin current. 

Let us consider a single vortex in a film in more detail. The Hamiltonian of local spin is 
modeled as 

S 2 d r 

H v = — d 2 x[J[{V6) 2 + sin 2 6{V(t>) 2 } + K± cos 2 9}. (91) 
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Since we cannot obtain the vortex solution analytically, we approximate it as 



1 V , 7T 
tan -1 - ± -, 
x 2' 



f r > A v 

-^(l-e-W) r < Av 



(92) 



where r = ^ x 2 + y 2 and A v is the size of the vortex core. Then the gauge field is given as 
Af ~ ±i x i(*±w) e -r 2 j 'A v 2 — neglecting the small contribution from outside the core. 
Choosing the current direction as x, we obtain the Fourier transform as (A is the area of the 
film) 

At(q) = ±m^(v^(l - 2q 2 X v 2 ) + 2iq\ v )e-^. (93) 

The resistivity due to the core is then calculated as 

vr 2 m 2 k F A 2 \ v 2 2C 2 kF 2 K i 
Pv 2 e 2 n 2 A 6 ' 

where we used the fact that resistivity is dominated by the contribution from q < k F + — k F - = 
2k F C- The nonadiabatic force is then given as 

F na ~ J -h ( {k F Kfde^ 2k ^\ (95) 
On the other hand, the force due to the topological Hall effect is given by 

F ad ~ J -h (96) 

for a vortex with n v = 1. 

Let us estimate the magnitude of the Hall effect. Hall resistivity is given as 

2ttS 2 Pti v 

The Hall conductivity is given by a xy = a o j+(^j^ ( a = Po 1 is the Boltzmann conduc- 
tivity) and hence the ratio of the Hall current to applied current is obtained as 

— = Pxy/pO- (98) 

Jo 

Thus, the deviation of the electric current becomes significant if the ratio ^ ~ 2irS 2 Pn v ^^ 
is of the order of unity. The deviation of current because of the Hall effect would be important 
in clean samples. 

In conclusion, we studied the nonadiabatic correction to the torque acting on local spin 
arising from the electric current, treating conduction electrons fully quantum mechanically. 
The nonadiabaticity we considered here is that due to the fast-varying spin texture. (The spin 
relaxation of conduction electrons results in the deviation of electron spin from adiabaticity, 
and the effect is sometimes also called nonadiabaticity. This effect was considered in our 
preceeding paper. 11 ) The nonadiabatic torque was shown to be nonlocal in space, to have an 
oscillation with a period of the order of 2k F , as was argued by Waintal and Viret 14 and Xia 
et al., 15 and to be of quantum origin. This oscillation torque is represented as a force acting 
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on the whole spin texture, and thus it represents the effect of momentum transfer. 1 ' 6 Note 
that this nonlocal torque is totally different from the local /3-term torque, which was originaly 
introduced to simulate momentum transfer force, 8 although they both contribute to a force. 
We found another type of force, topological or adiabatic force, which remains finite in the 
adiabatic limit and is proportional to the vortex number. These two forces were shown to be 
the counteraction of electron transport properties such as resistivity and the Hall effect. 

The authors are grateful to Y. Yamaguchi, T. Ono, M. Yamanouchi, H. Ohno, M. Klaui, 
Y. Nakatani, A. Thiaville, A. Brataas, R. Egger, M. Thorwart, J. Ieda, J. Inoue, S. Maekawa 
and H. Fukuyama for valuable discussion. G.T. is grateful to the Center of Advanced Study, 
Oslo, for its hospitality during his stay. 

Appendix A: Details of calculation of spin density 

We present, in this section, details of the calculation of spin density. Equation (25) is 
evaluated to the lowest order in Qq as 

Si 



Si + S 2 + S 3 + 5 4 + S 5 



iokq i 



B / 



ujkq 



S 3 

Sa 



u)kq 



ujkq 

[(£fe-2 )Ti J 2 5fc+2, ±A , - 9l-i, T:0J (9k+^±,J 2 ~ (tfjUf^J^L+f^u, + 9k-z,=F,u,(9k+%,± f 



u)kq 
u)kq 



\2/ a 



(9l 



? ~ (9l 



f(9] 



where we used partial integration with respect to ki and the identity 

m ^2d ki Bi = B . 

i 

The correction to the current vertex, eq. (27), is calculated as 



(A-l) 



(A-2) 



,5± 



om) (b) = E 



Lukq 



Bo 
fin 
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= -S 2 -S 3 -S, + Y. B ^'W k - h ^9l +i ,±^ 

u)kq 

Thus the current-driven part is given by 
^(x,t) = s ± (x,t)^+S ± (x,t)^ 

u)kq 

Lokq i 

X 
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(A-3) 



u)kq 



u)kq 



u)kq 



A.l Dominant term 



(A-4) 



We consider, in this paper, a clean limit, (e^r) 1 — > 0, where r is the elastic lifetime due to 
impurities, which we assume to be spin- independent. In this limit, spin density is dominated 
by Si. Thus 

S^\x,t) = £J>i» 

uikq i 



X [\ k+ 2Ji 5 fe-f'=F,wffe+2,±,a;ffc + 2 i±)a; + (& ~ ^ J . 5 fc-f , T ,a,5 fc-2 )TjW 5 fc+2 >±j , 



(A-5) 



Since ksT jep <C 1, we can replace /'(w) — — <5(w), and hence 



kq i 



X + f)^fe-2, T 5fe + 2, ± 5fe + 2 i± + (fc " |) .5fc-2, T 5fe-2 , T 5| +§ ,± 



where 



= 5fc-2 =p w =0' Using the identities 



9k,a9k,a = iT i9k,a ~ 9k,a) 



9l-i T 9l+i± 



e k+ i - e k _i 



T 2A + i 7 ( ^-f.T~^+f. ±J 



(A-6) 

(A-7) 
(A-8) 
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9 k - h ^ k+i ,± = - efc+ ,- 6fc _ aT 2A ( ^-f^-^ +f ,J> (A-9) 



where 7 = ^, we obtain, in the limit of r — > 00, 



2vr ^ ^ 2ml/ 

kq i li 



2 /y efc+g - £fc ± 2A ' 2 efc+ g — efc =F 2A 

-m Qi (gl^JMk + — £fc T 2A) - 5g >± J„(-2(fc + |))<5(e fe+9 - e fc ± 2A)) 

(A10) 

Appendix B: Estimation of S11 and sj_ 

We carry out the estimation of sy = ^ y^ ± e^^s^ and sj_ = ^ J2± ^fie^^s^. From eq. 
(A-10), we obtain 



1 



-ieEjT 



2 EEE 

± fcq i 



-iqx 



9k,± ~ gfc,± 
e fc+ 9 — £fc ± 2A 



-i7r ft <S(e fc+ , - e fc ± 2A)^ (fc + |) (^, ± AJ + Sfc^f) 



(B-l) 



where A±(a;,g) = e^^i^) and 



± fcq J/^ 



iqx 



m 



ki(k + 



Af _ A1 9 k , ± 9 k± 



+iirqi5(e k+q - e k ± 2 A)— (k 



m 



(B-2) 



From the symmetry argument, we see that kj (k + — > -j-fc ■ (fe + | ) and (ft (fc + |) 
■ (k + f). We thus obtain 



er 



£ e- iqX i(9 k ,± ~ 9%,±) \(E ■ A«)fc ■ (fc + f ) 



fcq± 

±7r(£ • A^g • (fe + |) <5(e fc+9 - e fe ± 2A) 



efc+q - e fc ± 2A 



kq± 



(B-AW)fc.(fc + |) 



£fc+ g - e fc ± 2A 
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A a »(x,q) EE 



e s4> (x) ■ A M (g) 
e^x) ■ A M (g). 
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(B-3) 



(B-4) 



Thus we obtain eq. (34). 
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